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Abstract. We report on a 30 ks observation of the high-mass 
X-ray binary pulsar OA01657^15 performed by BeppoSAX 
on September 1998. The wide band spectrum is well fit by 
both a cutoff power law, or a power law modified by a high 
energy cutoff, plus a fluorescence Iron Une at 6.5 keV. The two 
models are statistically equivalent. The inclusion of a cyclotron 
resonance feature at ~36 keV - coiTesponding to a magnetic 
field strength 3.2 (1 + z) x 10^^ G, where z is the gravitational 
redshift - improves significantly the for the cutoff model 
but only marginally for the power law plus high energy cutoff 
model. The statistical significance of our data is not adequate to 
discriminate between the two models, and even the Normalized 
Crab Ratio technique, successfully used to pinpoint cyclotron 
features in the spectra of other X-ray pulsars, is not conclusive 
in answering the question whether the feature is real or it is an 
artifact due to an improper modeling of the continuum used to 
fit the data. 
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1. Introduction 

X-ray binary pulsars (XRBs) display a common signature in 
their X-ray spectra: the presence of a high-energy cutoff that 
makes their spectra rapidly drop above ~40-50 keV. There are 
only few exceptions to this general behaviour, and OA01657- 
415 is one of them: indeed, this source displays one of the 
hardest X-ray spectra amongst XRBs, with no steepening up 



neutron star magnetic field. For OA01657^15 we therefore 
expect Bi2 10: only another XRB, A0535H-26, is known to 
possess such a strong magnetic field ( prove et al. 1995 ). 

OA01657^15 was discovered by the Copernicus satel- 
lite ( Polidan et al. 1978 ), and soon after its discovery a 38 s 
pulsation was detected by the A2 experiment aboard HEAO 1 
dWhite & Pravdo 1979| ). The optical counterpart of OA01657- 
415 was first identified with the spec troscopic binary V861 
Sco, but subsequent observations (Armstrong et al. 1980; 
Parmar et al. 1980| ) produced positional eiTor boxes that ex- 
cluded this star as the optical counterpart that is, up today, 
still unknown. Since 1991 OAO 1657-4 15 has been moni- 
tored continuously by the BATSE experiment aboard CGRO. 
This allowed the discovery of X-ray eclipses by the stellar 
companion and the determination of a 10.44 d orbital pe- 
riod ( [Chakrabarty et al. 1993 ). The distance to the source is 
not known, but a lower limit of 1 1 kpc has been determined 
by assuming that the neutron star is rotating at its equilib- 
rium period. Because of its low galactic latitude, the low en- 
ergy (<2 keV) spectrum of OA01657-415 is very absorbed, 
with Nh ^ 10^'^ cm~^. GINGA observations of the source re- 
vealed the presence of a soft excess, below 3 keV, modeled with 
a 0.3 keV blackbody ( |Kamata et al. 1990| ). The OA01657- 
415 X-ray spectrum also displays a fluorescence Iro n line at 
~6.6 keV dWhite & Pravdo 1979] [Kamataetal. 1990|). 



to 100 keV (Byrne etal. 1981). In those cases in which cy 



clotron line features are observed, the cutoff energy seems 
to be related to the cyclotron line energy Ecyc through a phe- 
nomenological relation found by Makishima & Mihara (1992): 

E,y, ~ (1.2 - 2.5) • E,. Because E,y, is related to the neu- 10 keV; [Boeflaet al. 1997bD , and the two mechanicafly colli 



2. X-ray observation and data analysis 

OAO 1657^ 15 was observed by the Narrow Field Instru- 
ments (NFIs) aboard the Italian-Dutch X-ray satellite Bep- 
poSAX ( jBoeflaetal. 1991?^ from 1998 September 4 15:41:46 
to September 5 07:20:35 UT. All the NFIs performed nomi- 
nally during the observation, namely the two imaging instru- 
ments LEGS (0. 1-10 keV; iParmar et al. 19971) and MEGS (1.5- 



mated instruments HPGSPG (4-180 keV; Manzo et al. 1997) 



tron star magnetic field strength B12 in units of 10 G by the 

relation E,y, = 11.6 B^2 (1 + z)'^ G, where z is the gravita- PDS (15-200 keV; [Frontera et al. 199% Data were col 



tional redshift, a measurement of Ec can give an estimate of the 
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lected in direct mode, providing information on time, energy, 
and position. Good data were selected using default criteria. 
The net exposure times for the four NFIs were 9.2, 26, 11.7, 
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Fig. 1. OA01657-415 count rate spectrum (plus signs) and best-fit continuum model (histogram), togetiier with the fit residuals. 
On the left the EC model was used, on the right the HEC model. The CRF normalization has been set to zero, in order to show 
its profile in the residual panel 



and 10.8 ks for the LEGS, MECS, HPGSPC, and PDS, respec- 
tively. These differences are due to the constraint that the LEGS 
must operate in satellite night time, and rocking of the colli- 
mated NFIs. For the imaging NFIs data were extracted from 4' 
circular regions centered on the source position. 

The background for the imaging NFIs was estimated from 
standard background blank field measurements. Because of its 
low galactic latitude, and therefore possible contamination by 
background sources, we also extracted spectra from two blank 
4' circular regions offset by 11' with respect to the source. Our 
result on the spectral analysis does not depend on which back- 
ground was used: for example, the 6-7 keV MECS count rate in 
the offset backgrounds is less than 1 % of the source count rate 
in the same band. For the non-imaging NFIs the background 
was monitored by rocking the collimators off-source by 3° for 
HPGSPC (one-side rocking), and 3? 5 for PDS (two-side rock- 
ing). The off-source fields were checked for the presence of 
contaminating sources: one of the PDS background fields was 
found to show an excess count rate with respect to the other, a 
clear signature of contaminating sources, and therefore was not 
used for the estimation of the source background. 

The contribution of the Galactic Ridge (GR) emission was 
estimated by using the results from RXTE observations of the 
galactic plane (Valinia & Marshall 1998). We found that the ra- 
tio between the GR and the source flux is always ^ 10^"^, 
therefore we exclude any contribution of the GR emission to 
our spectral data. 

The spectra of each BeppoSAX NFI were first fit separately 
to derive the best model able to describe the broad band spec- 
trum. We find that the LEGS and MECS spectra can be de- 
scribed in terms of an absorbed power law plus a Gaussian 



line in emission, while the PDS spectrum clearly presents a 
deviation from a pure power law, with a steepening above 
^-^20 keV. We therefore used a power law, modified at low en- 
ergy by photoelectrical absorption and at high energy by a cut- 
off, to describe the broad band spectrum. Variable normaliza- 
tions among the NFIs, found in agreement with those obtained 
from the fit to the Crab spectrum, were introduced to take into 
account known calibration uncertainties. Both an exponential 
cutoff (EC) and a high-energy cutoff (HEC; [White et al. l"983| ) 
yield a good fit, while a thermal bremsstrahlung did not fit our 
data. Also the smoother Fermi-Dirac cutoff (Fanaka 1986) did 



not adequately describe the high energy tail. An F-test shows 
that the probability of chance improvement (PCI) of the be- 
tween the EC and HEC models is 78%, therefore they are statis- 
tically equivalent. In Table |l| we present the best-fit parameters 
for both the two models. 

Concerning the low (^3 keV) energy part of the spectrum, 
our data do not sho w any soft excess as observed by GINGA 
( Kamataetal. 1990 ). The Iron line equivalent width is 360 ± 
280 eV, in agreement wi th the GINGA result bu t lower than the 
HEAO 1 measurement (White & Pravdo 1979). This confirms 



the possible contamination by the GR emission of the latter, as 
discussed by Kamata et al. (1990). The Iron line energy found 
by BeppoSAX is lower than that measured by GINGA (6.60 ± 
0.08 keV) and HEAO 7 (6.7 ± 0.1 keV). The unabsorbed fluxes 
in the bands 2-10 and 10-100 keV are 4.7 x 10"^", and 1.8 x 
10^^ erg cm^^ s^^, respectively. The total 0.1-100 keV X- 
ray luminosity is 2.5 x 10'^^ erg s~^ d^pc' where dkpc is the 
OA01657-415 distance in kpc. 

From the form of the residuals above 20 keV we were 
led to add an extra component to our spectral model, that 
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Table 1. Best fit parameters for the two models that fit the OA01657^15 continuum: a cutoff power law and a power law plus 
high energy cutoff. For both the models we show the best fit parameters without the cyclotron resonance feature (CRF) and with 
it. Because the CRF width is not well constrained by the fit, we present also the best-fit parameters for a CRF width fixed at 
lOkeV 



Parameter Value" 







Cutoff 




Power Law + Highcut 




w/o CRF 


w/CRF 


w/CRF 


w/o CRF 


w/CRF 


Nh (10^^ cm"^) 


12.1 ±0.4 


12.0 ±0.4 


12.1 ±0.4 


12 7+°* 

-0.2 


-19 7+0.3 
^^^•'-0.4 


a 


0.41 ± 0.04 




0.43 ±0.04 


0.83 ±0.04 


r, oq+0.03 


Ec (keV) 


16.4 ±0.5 




17.5 ±0.7 


13.0l«i 




Ef (keV) 








91 n+0.6 
^^•^-0.4 


99 9+0.8 
^^•^-0.7 


Epe (keV) 


6.48 ± 0.04 


6.48 ± 0.04 


6.48 ± 0.04 


6.47 ± 0.04 


6.47 ± 0.04 


(TFe (keV) 


r, qn + 0.08 
U.OO_g 


'-'■'^^-0.07 




0.42 ±0.05 


42+°-°'5 
'^•^^-o.oe 


EcRP (keV) 




30l25 


36±3 




44+11 


TCRP 




o.isin^ 


0.12 ±0.05 




11+006 


WcRF (keV) 






10 (fixed) 






Xl (dof) 


1.129 (572) 


1.088 (569) 


1.094(570) 


1.127 (571) 


1.117 (568) 



" Uncertainties at 90% confidence for a single parameter 



we interpret in ter ms of a cyclo tron resonance feature (CRF). 
Both a Lorenzian (Mihara 1995) and a Gaussian in absorption 
(Soong et al. 199C) improved the fit, and here we will present 
results obtained with the Lorenzian form of the CRF. The inclu- 
sion of the extra component have different significance for the 
two continuum models. Indeed, an F-test shows that the inclu- 
sion of the CRF improves significantly the for the EC model 
(the PCI is 2.8 x IQ-^), while the PCI for the HEC model is 
higher (1.1 x 10~^). In Table |l] we show the best-fit parame- 
ters of the CRF, and in Fig. [l| we show both the EC and HEC 
best-fit spectral models, with the residual CRF Lorenzian line 
profile shown in the residual panel, obtained by setting the line 
normalization to zero. We can see very well how the line is sig- 
nificant and very broad by using the EC model, and with much 
less significance with the HEC model. Because the line width 
is not well constrained by the fit, we also present in Table |l] the 
best-fit parameters obtained by fixing the CRF width at 10 keV. 

In order to check if the CRF is real or an artifact due to an 
improper modeling of the continuum, we calculated the ratio 
between the HPGSPC and PDS count rate spectra with the Crab 
pulsar count spectrum. This ratio has the advantage of mini- 
mizing calibration uncertainties. In Fig. ^ top panel, we can 
see that this ratio does not show any evident change of slope, 
as observed in other pulsars (see e.g. prlandini et al. 1998b ; 
Dal Fiume et al. 1998). By multiplying the Crab ratio by the 
functional form of the Crab spectrum (Fig. ^ second panel), 
i.e. a power law of index 2.1, and diving by the functional 
form of the continuum (Fig. ^, third and fourth panels), as 
obtained from our spectral fits, we are able to emphasize any 
deviations from the continuum. This Normalized Crab Ra- 
tio technique has been successfully used to pinpoint CRF s 
in the spectra of other X-ray pulsars ( Dal Fiume et al. 1999 ). 
The third panel of Fig. ^ shows a broad line, modeled with a 
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Fig. 2. OA01657-415 NormaHzed Crab Ratio. Top: Ratio be- 
tween the HPGSPC (black dots) and PDS (grey triangles) count 
rate spectra with the Crab pulsar one. Second Panel: Crab ra- 
tio multiplied by the functional form of the Crab spectrum, a 
featureless power law of spectral index 2.1. Third Panel: Crab 
ratio times E^^-^ divided by the EC continuum. Fourth Panel : 
The same as the third panel but divided by the HEC continuum 



Gaussian line centered at 43^5^ keV and dcyc = 28^g^ keV 
(xVdof ^ 17.7/30). 

The narrow feature at ^^60 keV present in the PDS Crab 



ratio has a 



keV if modeled with a Gaussian in 53- 
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70 keV (x^/dof = 0.18/4). This value is marginally com- 
parable with the PDS energy resolution of 8.6 keV (FWHM; 



Frontera et al. 1997) at 62j^2 keV. The inclusion of this feature 
in the broad-band fit yields a x^/dof of 618.19/567. An F-test 
shows that the PCI of the is 17.5%. For these reasons we do 
not consider it as a real feature. 



3. Discussion 

The main result obtained from the BeppoSAX observation of 
OA01657^15 is the clear detection of a steepening of the 
spectrum above 20 keV, and that this drop occurs in two steps. 
Indeed BeppoSAX observations of XRBs have shown that there 
are actually two changes of slope in the high (E<;10 keV) en- 
ergy part of their spectra: a first change of slope occurs in 
the ^10-20 keV range, while a second steepening occurs for 
higher energies. This is partly expressed in the form of the high- 
energy cutoff (White etal. 1983), that has been used to em- 
pirically describe the two-step spectral drop in these sources. 
In the XRB Vela X-1, for example, the first change of slope 
occurs at ^24 keV and it can be misinterpreted as a CRF 
dOrlandini et al. 1998a|). For OA01657-415 we think we ai-e 



observing the same phenomenon: the high energy part of the 
spectrum presents two changes of slope: if we model the con- 
tinuum with only one steepening, then the fit requires the pres- 
ence of a second one, that can be modeled with a CRF. On the 
other hand, if we model the continuum with a two-step change 
of slope (the HEC model) then an additional CRF is not re- 
quired to describe the spectral drop. The fact that we observe a 
broad CRF is due to the fact that the OA01657^15 spectrum 
is very hard, therefore the steepening is smooth. On the other 
hand, in the case of a softer spectrum, like in Vela X-1, the 
steepening is much more abrupt, and therefore tiie false CRF is 
narrower. 

Lacking a model of the continuum able to describe this be- 
haviour, the only way to discriminate a real CRF from an arti- 
fact due to an improper modeling of the continuum is through 
the Crab ratio, which is model independent. In Fig. || we show 
the Normalized Crab Ratios from PDS data obtained for a sam- 
ple of XRBs observed by BeppoSAX. A comparison with Fig. ^ 
clearly shows that the OA01657^15 Crab ratio does not pin- 
point the CRF as it does for all the other XRBs. The fact that 
Ecyc depends on the adopted continuum seems to rule against 
a real CRF, but it is intriguing that the Makishima & Mihara 
relation is satisfied for the OA01657^I5 CRF. Furthermore, 
the relation between Ecyc and its FWHM, clearly shown in 
Fig. 1^ where the CRF width increases as Ecyc increases, and 
explained in terms of Doppler broadening of the electrons re- 
sponsible for the resonance (Meszaros 1992) is fulfilled if we 
assume for the CRF width a value of 10 keV, which is the one 
obtained from the HEC fit, and can be considered as a lower 
limit for the width. In this case the OA01657^15 surface mag- 
netic field strength would be (3.1±0.2) x lO^^ (^i^^^ q Qn the 
other hands, this relation would not be satisfied by the broader 
CRF. Pulse phase spectroscopy should be able to resolve this 
ambiguity. However, longer observations may be needed. 
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Fig. 3. Normalized Crab Ratios from PDS data obtained 
for a sample of XR Bs observed by BeppoSAX (see e.g. 
Dal Fiume et al. 1999 ). The continuum drop is modeled with 
a HEC for all the sources but Cen X-3, for which an EC model 
is used 
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